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ABSTRACT 


The basic aim of this thesis has been to design 
regulator and stabilizer for asymptotic rejection of 
constant disturbances in a local load connected to a 
generation bus from the infinite bus and obtain improved 
dynaiiiic behaviour of a given power system. An appropriate 
linearized model of the given system is developed. Optimal 
modal control theory is used to optimize the controller 
parameters. The system performance with various contro- 
ller schemes are discussed. 



CHAPTER 1 


IHTRODUCTION 

Advances in the technolcgy of energy conversion, 
characterized hy increased sizes, pressures, and temperatures 
of steam generating units and large hydel units have resulted 
in bulk power generation and transmission. In recent decades 
substantial increase in reliability and reduction in cost are 
being achieved by interconnexion of these generating units. 

The fundamental problem in the design and operation of such 
systems are those associated with system stability, economic 
operation, reliability of supply to the consumers and security. 
Furthermore system voltages and frequency are to be constrained 
to be with-in specified limits. Many of these problems can 
be posed as problems in the area of automatic control. 

Modern control theory offers various techniques for 
the solution of such problems. Among these the application 
of Modal control theory T results in simplified and effici- 
ent techniques to develop designs of regulators for complex 
systems of large order. Simple design techniques are available 
in literature j 2 / fer realizing optimal modal controllers. 

The design of regulators is nearly always based on 
linearized system models [_ . Such models are always in - 

accurate. Inaccuracies may arise from auproxiraations made 
in the theoretical description of the system, from lineariz- 
ing a nonlinear model about an operating point, from simpli- 
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fying a higher order model to obtain one amenable to analysis 
from errors in parameter identification and from many other 
sources. In addition the systems to be controlled are sub- 
ject to disturbances, usually unmeasurable, and to changes in 
inputs and set points. 

Linearized models are generally used to design power 
system controllers. Often, reduced order models of these 
linearized models are used to simplify the computations in- 
volved. 

In practice power systems are subject to a large class 
of disturbances. Disturbances of im'^ulsive nature consist of 
those due to lightning strokes and instantaneous short cir- 
cuits in transmission lines. Switching on or off of loads, 
change in impedance due to removal of faulted lines and many 
sustained faults contribute to step ( that is, constant ) dis 
turbance. Beyond these the system encounters arbitrary dis- 
turbances due to hunting in the generators, arcing ground in 
transmission lines, power swings etc., many of which can be 
approximated in the class of polynomial disturbances for the 
sake of analysis and design. 

With such a system with frequent disturbances of a 
large class the aim of a satisfactory regulator design con- 
sists of determining controllers that yield good transient 
response and maintain output regulation, keeping in mind the 
simplicity and economy of the scheme. 
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G.D. Johnson |_ 4 f considered the case of external dis- 
turbances wliich are not directly measurable but can be re- 
presented as m^^- degree polynomials in time with unknown 
coefficients. The work has been extended by the author 5 { 

to acoamodate the case of unmeasurable vector disturbances 
ii li 

satisfying a. f - degree linear differential equation with 
known coefficients and vector control. Davison and 5rsith ’i_6j 
have found conditions for existance of a feed-back control 
system of minimum order for closed loop pnle placement and out- 
put regulation in presence of constant external disturbance. 

The same authors give design methods ^ 3 j linear regulators 

of both integral feed back and feed forward types for linear 
time invariant multivariable systems with constant disturbance. 
Davison has considered the case of assymptotic rejection 

of any measurable disturbance by use of feed forward control 
for multivariable linear time invarient systems. Davison and 
Smith j s'] have also obtained explicit feed forward and in- 
tegral feed back controllers fcr linear time invariant multi- 
variable plants, subject to both measurable and unmeasurable 
constant disturbances, so that the outputs are regulated to 
preassigned set points. 

This work considers a power system, consisting of a 
single generator oonneoted to infinite bus system with a local 
load at the generation bus. Controllers are designed for 
constant disturbance in the local load so that the infinite 
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bus system gets decoupled from the loc^l load disturbances in 
steady state and asymptotic regulation of voltege is obtained 
at the generation bus. Further more transient response of 
the system is improved by appropriate closed loop pole place- 


ment . 
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CHAPTEE 2 

PROBLEM FORMULATION 

2,1 INTRODUCTION 

The system under consideration ( Pig, 2.1 ) consists of 

a generator connected to a local load at the generation bus 
and a tie line connecting the generator to the infinite bus. 

The fluctuations in the local load get reflected on the genera- 
tor as well as the infinite bus system. The system, though 
considered to be infinite is in practice a finite one and the 
disturbance in the local load is transmitted even to the remote 
ends of the whole system, there by, bringing down the stability 
limit and the maximum power flow in any part of it . 

The problem comprizes of regulating the generation bus 
voltage and the power angle so that the disturbances in the 
local load are constrained to the local generator only there by 
improving the overall system performance. The class of cons- 
tant disturbances is considered for the controller design; 
impulsive disturbances are taken care of by the design of a 
stable system with good transient response. 

Optimal modal control theory is used for designing 
three types of controllers , namely 

(i) Peed forward, controller 

(ii) integral feed back controller 

(iii) Combined feed forward and integral feed back controller 
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The efficacy of the various schemes is studied. 

The above controller schemes are briefly discussed 

below. 

2.2 INTEGRAL PEED BACK AND PEED FORWARD CONTROLLERS \ 3 I , 

[4]: " 

Consider the system* 

• 

X = AX + BD + DW 

S; (2.1) 

Y = CX + EU + PW 

where X is the n-state vector, D the r-plant input 
vector, Y the m-output vector and W is the q-oonstant distur- 
bance vector. A, ..., P are real constant matrices of appro- 
priate dimensions. Assume without loss of generality that the 
system S is expressed in the minimal form, i.e. 

A B d"! 

Rank j=n+m (2,2) 

_C E P i 

It is desired to find a linear time invariant contro- 
ller for S such that Y(t)— ^0 as t — ^ and the closed loop 
system is stable. 

Necessary and sufficient conditions, for a solution to 
the problem to exist are? 

(i) (A, b) is stabilizable 

-rA ■■ B~\ . ' ■ 

(ii) Rank | | = n + m (2.5) 

■■ Lc ■■ Ej; ' 

It may be noted that the condition (ii) implies that 

r m.' , 
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We may further stipulate that the transient response 
of the closed loop system should he satisfactory. This will 
require assignment of closed loop eigen values to appropriate 
locations. This condition will require controllahility of the 
dominant eigen values of the open loop system. 

Assuming the above conditions to hold, let the desired 
final (steady state) values of X, Y and U be X, Y and TJ res- 
pectively. 

In steady state X = 0, and Y = n (by design require- 
ment). 

In steady state equns. (2,1 ) reduce to 


AX + BU + DW = 0 
CX + EU + P¥ = 0 


(2.4) 
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A o' 

r 

Let 

A 



and B i | 




C 0 1 
— ^ 

i- 

Prom 

e qun . 

(2 

• 6 ) we have 


P 

-- 

A P 

+ B Q 


Assuming a linear control law., 

Q = K P 

for the system described by equn. (2.7) we obtain 

V = K^Z + KgY (2.8) 

where 

K = jK^ I K^j 

Integrating eqn. ( 2 . 8 ) and substituting for V and Z . 
we obtain 


TI 


Prom eqn. 

TJ 


= K^X + Kg j Y(s)ds - K^X + U 
o 

(2.9) and (2.5) we obtain 
= K^X + Kg ^ T(s)da + K^W 


( 2 . 9 ) 


(2.10) 


where 


K !k 

? L 1 


t 

f ^ 

I 

f 


,T 




,T 


D 

P 



The state feed back K^X is used to place properly 
the dominant eigen values of the system matrix A to achieve 
desired transient response. 

The cl if Ye rent control schemes that can be implement- 
ed. are s , 
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(i) Feed forward control 

U = K,X + K^W (2.11) 

15 

(ii) Integral feed 'back control 

t 

U = K.X + K. j Y(s)ds (2.12) 

6 

(ill) Combined feed forward and Integral feed back control 
a = EX + K„ f Y(s)ds + K W (2.13) 

' o , 

The control schemes (i) and-(iii) necessiate that the 
disturbances be measurable. 

Optimal modal control theory can be successfully used 
to evaluate the feed back matrix. 

2.5 SYSTEM MO BEL 

Our system comprizes of a synchronous generator with 
an exciter and a governor and primemover. The generator is 
connected to the infinite bus system by a tie line with neg- 
ligible resistance. At the generation bus a load is connect- 
ed which is subject to disturbances.' The load is represented 
by constant impedance for the design and the dist'urlifences are 
assumed to be constant. 

A model linearized about an operating point is Obtain- 
ed for the above system. 

2.5.1; LohdRepresenta'fcion; 

Represent ihg the terminal net'wdrk by its Thevehin’s 
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where 


r'' + jx'' = EX + .iAX 


V' = AY + ,iBYo 


(2.14) 


y"* represents the equivalent voltage source 
1 1 

R + jX is the equivalent series impedance 
R^ + X(X + X,.) 

A A. L_ 

R'^ + (X + X, 


R^ + (X + 


2,5,2 Machine Equations \^\ J 

mmm m» mm mm mm mm mm is mm w^mm mm 


"V-fd “ ^ffd ^fd \d ^hd " \d ^d 



= X , i», 
md fd 

\d ^kd 

^d "d 


= X , i^, 
md fd 

^kkd ^kd 

- \d ^d 

T 
' <1 

^mq ^kq 

, X i 
q q 


^q 

“ ^kkq ^kq 

- V \ 


ij' 

■fd 

= — p V 
w,. • fa 

^fd ^fd 


®d 

= ” p y- A 

Wo ' a 

W 

- r X. - — 

a Wr: 

■ V 

^ q 


(2.15) 


O — -i— rs V' 4- 3? i 
■ w^ ^ ^ ^kd kd 
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Cos Cos ( o" + tan 'Sifk') 

= V Vo 

6"^ ■= o' + tan"'' B/A 


Taking the machine currents, load angle ' '{ ' arid mach- 
ine speed ' V.' ’ as the state variahles and eleminating the non- 
state variables from the above sets of equations and with 
appropriate substitution for the terminal voltages we obtain 
the equations: 
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^ffd ^ ^fd ^md ^ ^kd " ^md ^ ^d ~ '^o^^fd " ^fd ^fd^ 


X . P + X , p i, , - (X, + X^ )p i, = vr\ - \! + 

md fd md ^ kd '' d <3 OL ’ 

5< Y Sin (<r+ tan""^ B/A) + (r'' + r) i^ + 


B‘ 


— (X i, -(x'' + X ) i )! 
^ mt kq ^ q' q j 


X j p ix.j + X, , j p i, T - X - p i, = -w r, , i-i-j . (2.17) 
md ^ fd kkd ^ kd rad ^ d o kd kd 


X pi, - (X + x”*) p i = wj + 

mq ^ kq q <1 Oi. ' 


2 

B V 


Cos ( i + tan”"* B/A) + (r"* + r) i_ - 


— (X , i_, + X , i, , - (X, + x"”) i,)l 
w ^ md f d md kd d a ^ 


\kq P ^kq " \q ^ ^q "^o ^kq ^kq 


p /> = w - w 


P ^ “ ^i “ ^d + ^^rad: ^fd ■*■ ^md ^kd 


^d ^d^ ^q ^kq ” ^<1 ^d 

2.5.5 Voltage regulator |_ 9 I (fig* 2 . 4 ) = 

The voltage regulator is represented hy a linear 
model with an internal stahilizei* loop; 
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T"e 5 ®fa =■ ^<^,1 - ' 


_ V ) » 1 

•fd 


P ^ 
Y 


= k p E_ - V 
8 ^ id s 


( 2 . 18 ) 


fd 


•fd 

^md 


E 


fd 


2 , 5.4 Primemover and governor ! 11 j (fig. 2,5) J 

The primemover and the governor are represented hy a 

single time constant. The governing equation is taken to be 

^T. 


p AT. 

mg 


1 


A m 

rg • ing 


( 2 . 19 ) 


g 


2.3«5 System state equations! 


The above equations ( 2 .I 7 ) - ( 2 . 1 9 ) representing the 
system are a set of non-linear equations. These equations 
are linearized about an operating point. The operating point 
quantities are denoted by the subscript *0’ and the incremen- 
tal quantities about the operating point by the pre f'iz * 6 ' , 
Expanding in Taylor series and taking the first order approxi- 
mation the set of equations reduce to : 

g; X = 4.x + B.tr + D.w : ^ ( 2 . 20 ) 

w 1 1 1 1 

o 

where 

X is the statevector ; |a ^fd » ^kd ’ ■ - ^d * ^ ^kq J 

Iq , Ac3-' , AW , A T^ , 6 ; AY 

U is the control vector vI'aE: j AT^^^ | 

¥ is the disturbance vector ; [AR , AX, | 



^ffd 

X , 
md 

-X . 
md 

Sd 

1 

Sd 

-(x,+xS 

V d 0 

!X 

i md 

^kkd 

-X . 
md 


X 

-X +X 

mq 

q 

x, , 

-X 

kkq 

mq 




(1 0 


w T 
o ■ m 


w r 

o g 


w T 
o e 




T 7 


0 I 


The non-zero entries of and D. are : 


1^(1 , 1 ) = 

» 9) = ^fd/^md 

A^(2 , 3) = Rq + r 

AA2 , 4) = X 

^,(2 , 5) - -(3tJ+3C^) 

A ^(2 , 6 ) . -T^ 0 ,„: 

Al(2 , 7) = -(V’t;) ijoAo 

1^(5 , 2 ) . 

^<4 , 1 );- -x^ ' .V 


4/5 , 4) = ■ 

A^(6 ,7) =1.0 

^1^7 ; 1) = So 

^1(7 » " Sd ^qo 

A (7 . 5 ) = (X -X, ) i 
I'"' q 3 qo 

A^(7 f 4 ) = "Sq ^do 

S(Sv 5) = x^a Sdo-^<S-S^ ^ 

■ 4 .^( 7 :,: 8 .),.= 'i.o.:' 
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(4 , 

2 ) 

= -\d 



(8 , 

8 ) = 

= -1.0 


(4 , 

5) 

= X + x"* 
d 0 


"1 

(9 , 

9) = 

= -1 .0 

"1 

(4 , 

5) 

= r''+ r 

0 



(9 , 

10) 

= -K 

e 


(4 » 

6) 

= -^0^20 



(10 , 

- 9) 

= /r 

s 


(4 , 

7) 


+(x^+x;) 

^do 

)/Wo 

A. 

,(10 , 10) 


(9 , 





F ] 

c 


®1' 

1 ) 

= -^e ^ 

(10 , 1) 

= 

e 

s . 

B^(fl ,2) 


e 


D ^(2 , 1 ) = X 


I),(3 , 1) - Rso^^'^do-^ V20; V* V20^ 

I)lC2 . 2) = X2oC-V^il,0- X 1^„) - X^„(X 

D,(3 , 2) . X2„(X - =^10'* 1,5-^ ^o®20' 


where » 


Sin( £ +tan*'*B /A )A 


‘10 


o' o' o Ah2 , ,2 


\/^ 

V O 


+ B" 




-WAo + - 

M 

O 


B 


Cos^(tan”"*B /A ) 
^ o o' 


B 


10 


B Sin(<r +tan”'*B /A j sJa^+ B^ 
O O 0 o ^ o o 

"7 ^ “' 


A^ + 


^ : OAs^XiBn— B^/A 


C o s ( ^ + 1 a h ~ B / A ) 
o o o 


C <M 
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^ CosC £ +tan"^B /A ) 


'A^ + B 
o o 


r 2 2 

o o 


Cos(SL+tan”^B /A ) + 


B \/ A + B o _•! , \ 

.0^-2 a_cos^(tan B /A ) . 


Sin( ^+tan""*B^/Ao) 


fA^ + 

O o 


Cos(^^+tan* ’B^/A^) - 


fA^ + B^ P X 

L ° .-^ Oos^ftan ^B^/A^), 

xl 

o 


Sin(^ +‘ban"'' B^/A^) 


Coo = IM + ®o 
20 1/0 o 


2 + Sin(^ +tan"''B /A ) 

/-N O U w V 


B 2 b; 
0 0 

R ■ X 
o 


2(X + X^) B^ 


R X 
o 


2 

-2. (i _ A ) 

X ^ o-^ 


^9-- ( (1 :+ A ) (l + + -^o^ 


^20 = RT V- - X ' 0 


2.5.6 System_output _eQ-uations_i 


In the 


problem posea; the variehlee to he regulated are 

_ .ea the: load angle hetween;: the geaeratioh 


generation bus voltage and the 
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and the infinite hus voltages, 

[A<^2 

where | 

^2- ^ - 


Thus the output vector Y = 


'm 


A'Co = A& - 


in 


( 2 . 21 ) 


A-fit ^®n 

^ d e. q 

to to ^ 

Substituting for e^ , from eqns (2.16) , 

^ = Cos "'e /e. (fig. 2, 5) and expressing eqns. (2.21 ) in 

terms of the state, control and disturbance variables we 
obtain ; 


Y = CX + EU + PW 


( 2 . 22 '! 


where , 


C(1 , 1) = 


A(5 , 1) 

w E.- 
o 1 0 


, C(1 , 2) = 


V’^5 , 2) 

W E._ : * 

o 10 


C(1 , 5) = 


rE.-tan 7- + 

10 0 


'10 


x' A(h , 3) 

o • * 

w E. ^ 

0 10 


xli(5 , 4) , , o 




1 


c(l , 4 ) = 


x: ^'(5 , 5) 
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Vao < 

C(1 , 6) = 1.0 - g + - » 


10 


%®10 


0C1 • ’’'> - i®E,„ '-do 


(i. + *(5 , ^)'l 


o"10 
x! A(5 , 8) 


A(5 , 9) 


0 

td 

0 


x^ 

0 

A(5 , 

10) 


^o®10 

x^ 

0 

B(5 , 

1) 


w p. 

0 1 0 


x^ 

0 

D(5 , 

1) 

" 

0 

0 



, 0(1 , 9) = 


0^10 


, E(1 , 2) = -^ 


Bf'j , 2) 


%®10 


- .^,-(1,2) = 


x: D(5 » 2 ) X^Q 


r1 

'o 


E 


■w.E 


10 


10 


E 


10 


0(2 , 1) - C,o i(3 , 1) + =40 *'5 . 1) 

0(2 , 2) - Cjo i(3 . 2) + O 40 > 5) 


e, r'' 
do 0 


xl e 


0(2. , 3) = e.. "■ "30 


to 


+ 0 A(5 , 5) + + ^^40 

30 to 


C(2 , 4 ) = C 30 » 4 ) + a(5 , 4 ) 




G(2 , 6) 


e, las. v Go„+ ^/tn 

- V10 “30 -to “ “ 

■■ .to ■ . ■ ' 

0(2 , 7 ) . 030 (^ 0 .7); - 
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C(2 , 8) = C^Q A(3 , 8) + A(5 , 8) 

C(2 , 9) = A(5 , 9) + C^Q A(5 , 9) 


C(2 , 

10) 

■ =50 

A(3 , 

, 10) 

+ 0^0 A(5 

, 10) 


E(2 , 

1) = 

= 30 

B(3 , 

1) + 

°40 

b(5 , 

1) 


CM 

, 2) = 

‘ ^30 

B('3 , 

2) + 

°40 

B(5 , 

2) 


F(2 

, 1) = 

lis. 

®to 

*‘30 

°30 

D(3 

, 1) 


B(5 » 1) 

F(2 

, 2) 

lis. 

®to 

^30 

^30 

D(3 

, 2) 

+ 5^ ^40*”^ ’ 

2^°40 
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2.3 ♦? Numerical values ; 

The final siaie equaiions are 

X = AX + BTT + DW 
Y = CX + EIJ + FW 


( 2 . 23 ) 


where , 


A = w G7^ a 

0 1 1 

B = w GT'' B. 

o 1 1 

D = w G'"' 

o 1 1 


Machine parameters {. 9 ] : 


md 

as 1^0 

X = 0.6 

mq 

"d - 

1.2 X = 0.8 

G 

"ffd 

= 1.1 

x, , ^ 

kkd 

= 1.1 

X,, = 0.8 

kkq 

^dl“ 

0.2 Xj^= 0.1 

^kdl 

= n.i 

kql 

= 0.2 

Xf^=o.0011 

r = 

0.01 rj^^=0.02 

^kq 

= 0.04 

T = 
m 

6 sec 

k^ = 1.0 

w ~ 
0 

514.0 



f regulator 

parameters 

[ 9 ] 

• 

p . . 



T 

s 

= 2.0 

% = 0-5 

Ic = 
e 

40.0 = 0 

-04 

- ■ 


Pr’imeniover and g-overnor parameters 


X, = 0-5 

S 


Operating point 2. ' 

The operating point values chosen for the system are: 
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, Real power = 0.8 p.u. 

Reactive power = 0.6 p.u, lagging 
System frequency 'w^' = 514 .O rad /sec . 
Infinite "bus voltage = 1.0 p.u. 

Local load resistance =1.6 p.u. 

Local load reactance = 1.2 p.u. (inductive) 

The line parameters C?! s 


The tie line resistance is neglected and the reactance is 
taken to be O .25 p.u. 


With the above data the system matrices are found as 


9 
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CHAPTEF - 3 
DESIGN OP CONTROLLERS 

The discussions to follow present a systematic pro- 
cedure for the design problem. Firstly a check for the exis- 
tence of controllers of the kind discussed ia chapter 2.2 is 
made. This is done by checking whether the conditions given 
by eqns. (2,5) are satisfied, 

3,1 JJHECK FOR CONTROLLABILITY ; 

Conditions given by eqns. ( 2 , 3 ) require the system 
to be stabiliaable . The stabilizability is assured if the 
system is controllable. Thus we check for the controllability 
of the system under consideration. 

The controllability for a system with distinct eigeii 
r ' T 

values is obtained ! 1 \ iff no entire row of V B is identi- 

cally zero ; where V is the mstriy with the eigen vectors of 
T 

A as the columns. 

The required eigen vectors for the given system matrix 

r "i T 

A are calculated ; 12 j and the product T B found. 

Eigen values of the system matrix A : 

-19.85768566 + 0 313.15682556 

-19.85768366 - 0 313-13682556 

-86,76879025 
-27,88085959 
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-5-82955279 ‘ 

- 1.00000119 

-0.96654579 

-0.05847392 + j 0.52290527 
-0.05847392 - j 0.52290327 

-0.26044716 

I -0.00011917 + j 0.000821^67 O.OOOOOO44 + 3 0.00000461 

I 

-0.00011917 _ j 0.''00820<^7 O.OOOOOO44 - j 0. 00000^6 
-0.7^100«161 -0.00OO0484 

0.053248 0.0 

„ - 0.00011298 0.99999997 

V B = 

0.45947005 0.11687776 

0.44152659 + 3 0.52290349 -0.01207782 - 3 0.0653987 

0.44152659 - 3 0.52290349 -0.01207782 + 3 0.0655987 
-0.2394897 0.0 

Clearly the system is controllable. 

3.2 CHECK FOR THE RANK ; 

/ A b\ 

The matrix i for the system under consideration 

V C E / 

being a square matrix, its eigen values were calculated ['I2] 

It was found that non«. of the eigen values are zero. This- 
shows the rank of the matrix to be equal to the order 
(= n+m = n+r) of the matrix. 





- 19.857664 + j 315.13678 

-19.857664 - j 513-13678 

-86,768625 

-27 . 88145 
-3.3094957 

O.79I84762 + j 0.61496324 

0.79184762 - j 0.61496324 

-1.5901205 + j 0.54674523 

- 1 . 5'^01205 - j 0 .'^ 467''523 

-0,48^4694 

-O.3958B66I + j 0 . 4087921 ’; 

-0.30388661 - 0 O.4O879213 

Thus it is ascertained that controllers of the type 
discussed in chapter 2.2 exist for the given system. 

3.3 -PEED BIGK MATRIX CALCDLATIOR : 

The feed hack matrix is obtained by optimal modal con- 
trol theory as a sum of dyadic matrices. The dominant eigen 
values are devided into groups and for the pole assignment 
problem for each of the groups, dyadic feed back structure 
is assumed ';_ 2 ] . 

The objective function for minimization for any par- 
ticular group is taken ass 
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m 


J = ^21 

i=1 ■ 






(^. 1 ) 


where , 


T 

F = Os g is the feed hack matrix corresponding to the 


group . 


To shift the eigen values in group q after assignment of 

eigen values in the first (q-l) groups, the feed back aatrix 
is |_ 2 J : 


X«l) _ ^(<l) „(q)^ 




S 


* (q“i5»»««,p) 


(3.2) 


.where 


p represents number of groups 

3 (q) (q) 

i^i 


s 


(q) 


K. V. 
1 -1 


(3.3) 


1 


1 

Q 

"IT ( 

0 = 1 


j(q) 

' j 


>(q)) / xtih) 

^ ' 0=1 ^ ^ 


(3.4) 


(i = 1 ,...., 1 ) are reciprocal eigen vectors of 
■•1 q 

associated with eigen values (i = 1 ,...,, 1 ) and 

1 q 


4(q-i) 


i + B y g'i’ 

1=1 “ . 


.^T 

S' » 


(3.5) 


A is the system matrix and B the control matrix. 
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and (i = 

pending to the group o with 1 

Q. 

The (i = 1 ,... 

mode controllability condition 


) are the eigen values corres- 

q 

modes. 

m) are so selected that the 


B do, (i = 1,...., 1^) (3.6) 

is satisfied. 

Negative gradient technique ^21 is used for optimiza- 
tion and the resultant feed back matrix is calculated as ; 

(i) ■ 

K = X.. (3.7) 

i = 1 

3.3.1 Grouping criterion [_2j : 

To obtain low feed back gains, each pair of complex 

conjugate eigen values are shifted in a group. The real eigen 

values are grouped such that the vectors 7^*^^ fi = 1,....,l ) 

* 1 q 

are confined within as narrow a cone as possible. Further - 
more, the projections of the vectors (i-1,.,.,,1 ) on 

“1 q 

the space spanned by the columns of B matrix are as large as 
possible. 

3.3.2 Algorithm for optimal modal controller sign_using 


grouping criterion t 


(i) Determine the augmented system and control matrices 
A and B respectively of eqn. (2,7). 
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(ii) Determine the dominant eigen Values of the system 
matrix A and select new locations for them. 

(iii) If there are repeated eigen values displace them 
to different locations by applying. a suitable state 
feed backj so as not to e'f'fect much the mondominant 
eigen vslues. 

(iv) Set q = 1 . 

(v) Calculate the normalized left eigen vectors and 
eigen values of A 

"t h 

(vi) Form the q rroup using the grouping criterion. 

(vii) Initialize by checking the mode controll- 

ability condition. 

(viii) Set the objective function value J to a large 
value . 

(ix) Calculate ( 3 = 1 ,. *..,1 ) and g^^^ from eqns. 

(3.4) and (5.5). 

(x) Calculate the new objective function value. 

(xi) If 1 J - J , ,1 ^ (a small specified positive 

number) go to step (xiii) otherwise renew the objective 
function value. 

(xii) Change as + c( where 

c is step length and is the negative of the gra- 
dient of with respect to 2 j . Check 

whethef the mode controllability criterion for the 
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modes of the group is satisfied by the new 

vector. If it is not satisfied then reduce c apnropri 

ately to satisfy the criterion and go to step (ix). 

(xiii) Calculate the modal ^eontroller matrix and 

■ (xiv) If = p go to (xv); else increment q by 1 and 
go to (v). 

(xv) Calculate the modal controller matrix K from 
eqn. (?.?). 

(xvi) Stop. 

3 . 4 «'^ ® augmented system matrix : 

-19.857645 + 0 515.13678 

-19.857645 - j 515.15678 

-86.76852 
-27.88O85I 

- 3.839552 

-1 .0000009 
-0.96722747 

-0.0584884 + 3 0.52292036 

- 0.0584884 - 0 0.52292036 

0.0 
0.0 

3.4,2 New locations for the dominant eigen values s 

- 2,3 + .1 0.5229, - 2.5 - 3 0.5229, - 2 . 2 , - 1 . 9 , - 1 . 8 , 





Optimal 

mod^l controller matrix : 

. 1.5265 

1.4249 

-1.491 

5.615 - 

-1.5265 

-1 .4249 

1.491 

-5.615 

167.4762 

9.4585 

0.2088 

-9.5258 

-167,4762 

-9.4585 

-0.2088 

9.5258 


^•5 -22lroller_schemes_for_given_s7stem (2.25) : 

The matrices are computed from equations 

2.8 and 2*10. 


^1 = 


K2 = 


1.5265 1.4249 -1.491 5.615 -2.6862 
-1.5265 -1.4249 1.491 -5.615 2.6862 

180.9502 167.4762 9.4585 0.2088 -9.5258 

-180. 9502 -167,4762 -9.4585 -0.2088 9.5258 

66.0895 -22,7707 

-65. 0895 25.7707 


■0.9476 -22.9095 

0.9691 25.4795 
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CHAPTER - 4 

SYSTEM PEEPORMilNCE 

Constani: disturbance of 0.0625 p u in load resistance 
and. 0.0833 p u in load reactance is applied. The performance 
of the uncontrolled system and the system with various controllers 
are obtained by digital simulation. Runge-Yutta 4 "^^ order [l ll 
method is used for the solution of the problem. The system res- 
ponses for different cases are given in Rigs. 4.1»-..»4*5» 

The regulator designs reported are preliminary designs. 

The eigen value locations are chosen rather arbitrarily. The 
eigen value locations have to Tie so chosen that the regulator 
response is satisfactory and ho component of the regulator goes 
in to saturation for the usual system disturbances. This require- 
ment imposes the extent to which eigen values may be shifted in 
the left half plane. It is also desirable to limit the feed back 
gains. The responses of the system with the regulators seem to be 
a little slow. The performances can be improved further by shift- 
ing the eigen values further in to the left half plane. For a 
practical design satisfactory locations, for the eigen values can 
be chosen by trial and error considering the above-mentioned 


limitations 



CHAPTER - 5 


CONCLUSION 

A power system network consisting of a generator with 
a local load connected to an infinite bus system by a tie line 
was considered for regulator and stabilizer design. A stabili- 
zer was designed using state feed back to improve the dynamic 
characteristic of the system. Three types of regulators, (i) 
integral feed back , (ii) feed forward and (iii) feed forward 
in combination with integral feed back were designed for asym- 
ptotic rejection of constant disturbances in the local load from 
the infinite bus system and to obtain voltage regulation at the 
generation bus. An appropriate linearized model of the given 
system was developed and optimal modal control theory with group- 
ing criterion was used for the design of controllers with opti- 
mum controller parameters. The system performances were obtained 
with various controllers for a constant disturbance. 

In our design we have considered complete st«te feed 
back which needs all the system states to be measurable. 
Alternately we need to design an observer to obtain the un- 
measurable states. It may be possible to obtain a simpler con- 
troller by using incomplete state feed back 1^151 • 

The class of disturbances considered for the design 
consists of constant disturbances only. In an actual system 
the disturbances are of arbitrary nature. Design techniques 
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are available in literature for rejection of certain classes of 
disturbances* For rejection of polynomial disturbances the con- 
trol scheme becomes more complex as compared to that* for constant 
disturbances and the complexity increases with increase in the 
order of the polynomial. In case of a power system problem it 
may not be of much use to implement regulators for rejection of 
polynomial disturbances- 

The work hers considers a single generator to infinite 
bus through a tie line- The design technique can be extended 
for a more realistic network with mult imachine systems - Regulators 
can be designed with the guide-line given in the work for sharing 
a given load in various desired proport ions . and steady state load 
flow in different links of a power system network can be controll- 
ed accordingly- 

The representation of load considered here would be in- 
adequate in practical situations- The work can be extended by 
considering a more realistic load models. 

Wide range modal controllers can be designed r2l for 
reduction of sensitivity of controllers to change in operating 
point of the system- Accessible state feed back schemes can be 
tried to obtain simpler controller schemes. 

In the design the system performance equations are lin- 
earized about an operating point and controller schemes are de- 
signed to improve the performance of this system for small 
disturbances. It is worth studying the performance of nonlinear 
model with the linear regulator for large disturbances* 
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Fig. 2.1. Power System 
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Pig. 2.2. Thevenin equivalent 
circuit of terminal network 



Pig. 2.4. Voltage regulator representation 
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Pig. 2.5« Governor and primemover representation 
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4,2 SYSTEM RESPONSE WITH STATE FEEDBACIT 
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UNDO NT ROLLED SYSTEM 
RESPONSE (tr = 0) 


Scale 
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i Pig. 4.3(a) SYSTEM RESPONSE 
.• WITH STATE AND INTEGRAL PEED- 
"! BACK CONTROL 
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■i Fig. 4.5(b) SYSTEM RESPONSE WITH 
/ STATE AND INTEGRAL PEEDBACE: CONTROL 
' .t 

(H = K X + N j Y(s)ds) 
o 
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Pig. 4.4. (b) SYSTEM RESPONSE WITH STATE FEEDBACK 
AND FEEDFORWARD CONTROL 
(U « K^X + KjW) 
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Fig. 4.5(a) SYSTEM RESPONSE WITH STATE FEEDBACK, 

.FEEDFORWARD AND INTEGRAL FEEDBACK CONTROL 

(N - K^X + Kg ^Y(s)ds + K^W) 
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Pig. 4.5(0) SYSTEM RESPONSE WITH STATE FEEDBACK, 

FEEDFORWARD AND INTEGRAL FEEDBACK CONTROL 

(tl - K^X + Kg J Y^s)as + K^W) 
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